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We study generalized parton distributions in the impact parameter represen- 
tation, including the case of nonzero skewness ^. Using Lorentz invariance, 
and expressing parton distributions in terms of impact parameter dependent 
53^1 wave functions, we investigate in which way they simultaneously describe 

\ longitudinal and transverse structure of a fast moving hadron. We com- 

pare this information with the one in elastic form factors, in ordinary and in 
^ ■ /cT-dependent parton distributions. 



1 Introduction 

In recent years, generalized parton distributions (GPDs) PI El Ej have become a 
subject of considerable theoretical and experimental activity. Much of the interest 
in these quantities has been triggered by their potential to help unravel the spin 
structure of the nucleon, as they contain information not only on the helicity carried 
by partons, but also on their orbital angular momentum [2]. This is, however, only 
one of the aspects where GPDs can provide qualitatively new insights on how a 
hadron looks like at the level of quarks and gluons. A look at the variables on which 
GPDs depend immediately reveals that they simultaneously carry information on 
both the longitudinal and the transverse distribution of partons in a fast moving 
hadron. The physical picture becomes particularly intuitive after a Fourier transform 
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from transverse momentum transfer to impact parameter |4. Analogies with other 
imaging techniques have been highhghted in [5J: the Fourier transform also occurs in 
geometrical optics, with light rays corresponding to momentum space and an image 
plane to transverse position, or in X-ray diffraction of crystals, where the diffraction 
pattern has to be Fourier transformed to recover spacial information. 

Burkardt has shown that for vanishing transfer of longitudinal momentum (i.e. 
^ = in standard notation) GPDs transformed to impact parameter space have 
the interpretation of a density of partons with longitudinal momentum fraction x 
and transverse distance b from the proton's center It is natural to ask how this 
situation looks like at nonzero ^, which is relevant for most processes where GPDs 
can be accessed. This is the purpose of the present work. 

In section |2] we will see how the interplay between longitudinal and transverse 
degrees of freedom is dictated by Lorentz symmetry. We shall then study an explicit 
realization of this interplay in section IHl using the technique of [0] we will write GPDs 
in terms of the light-cone wave functions for the target, in the mixed representation 
of longitudinal momentum and transverse distance of the partons. We discuss some 
special cases of this representation in section and the question of resolution scales 
in section [S| Apart from GPDs, transverse information on partons in a hadron is 
also contained in parton distributions that depend on the tranverse momentum of 
the struck parton. We will compare the different nature of the information contained 
in these quantities in section El before summarizing our main results in section [3 
For definiteness we will present our discussion for quark GPDs in a nucleon here; its 
generalization to gluons and targets with different spin is straightforward. 



2 Lorentz invariance and its consequence 

2.1 From momentum transfer to impact parameter 

To describe the kinematics of GPDs we will represent a four- vector v in terms of light- 
cone coordinates = (v^ ± v^)/\^ and its transverse components v = {v^,v'^). To 
avoid proliferation of subscripts like "_L" we will reserve bold-face for two-dimensional 
transverse vectors throughout. The momenta and helicities of the initial (final) nu- 
cleon are p {p') and A (A'). We use Ji's variables x and ^ for plus-momentum 
fractions relative to the average momentum P = ^{p + p'), in particular we have 
{p — p')^ = C, {p + p')^- In the light-cone gauge = 0, the matrix elements defining 
GPDs are 



HvA = ^=l=/^e-^"^-(p',A'|g(0,-i^-,0)7+g(0,K'0)b'^) (1) 



1 



uip',\' 



H{x,^,t) 7+ + E(x,e,t) 



ia+°'{p' -p)a 
2m 



u{p, A) 
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for the sum over quark helicities, and 



HvA = -j=L= J^e^-''^^-{p',X'\qiO,-'^z-,0)^^^,q{0,lz-,0)\p,X) (2) 

(p' -p)+75 
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H{x,^,t)^-'^, + E{x,^,t) 



2m 



u{p, A) 



for the hehcity difference, where m denotes the proton mass and position arguments 
for fields are given in the form q{z) = q{z~^, z~ ,z). Note that these definitions do not 
refer to any particular choice of transverse momenta p and p' : Lorentz invariance 
ensures that H, E, H, E only depend on x, ^ and the invariant momentum transfer 



-t = ~{p'-pr = ^, + {i-e)^' (3) 

with 



D = ^. (4) 

That t depends on p and p' only through their combination in (0} is a consequence 
of its invariance under transverse boosts [ZIIH]. These are Lorentz transformations 
which transform a four- vector k according to 

k+ -^k+, k k - fc+v, (5) 

and k~ such that k"^ stays invariant. Here v is a transverse vector parameterizing 
the transformation. Since p'~^ = (1 — ^)P^ and p^ = {1 + O-^^j both vectors on 
the right-hand side of (jU are shifted by the same amount P+v, and their difference 
remains the same. 

In order to have simple properties of the matrix elements (0) and (j21) under these 
transformations, it is appropriate to choose the polarization states of the protons 
such that they have definite light-cone helicities [7||H]. The state vector and spinor 
for a proton with momentum components (p"*", p) is then obtained from the one with 
(p"*", 0) by a transverse boost ^} The matrix elements Hyx and Hyx are hence 
invariant under transverse boosts, and depend on the proton momenta only through 
the variables ^ and D, as do the GPDs H, E, H, E themselves. With the phase 
convention for proton spinors given in jU] we have 

^++ = n^- = H-^^E, 



^In contrast, the corresponding state with definite (usual) heUcity is obtained from the one with 
{p'^, 0) by a spatial rotation. 
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and 

7i++ = -H — = H - ^_ ^2 ^ ' 

~ ~ -\- ~ 

= (^+-)* = (7) 

We can now transform these matrix elements to impact parameter space, 

1 /""^ 



P(D^)Jo(|D||b|)(iJ-^£;), 



An 

1 6^ - ife^ /-"o 



-P(DV.(|D||b|)ME, (8) 



47r |b| 

with analogous relations for the other matrix elements. Notice that E and E parame- 
terize both proton helicity conserving and nonconserving transitions and thus appear 
twice, weighted by different Bessel functions Jq and Ji. 

We remark that there is some freedom of choice in introducing impact parameter 
GPDs. Instead of defining the Fourier transform with respect to D, one could for 
instance use the vector (1 — D (which coincides with p' — p in a frame where 
p' = — p) or the vector (1—0 D (which is equal to p' when p = 0). The corresponding 
GPDs are related to those in (jHI) by an overall factor and a rescaling of the impact 
parameter. Just as the different ways to parameterize the plus-momentum fractions 
of the partons [21 Ej, such a change of variables does not change the physics content 
of the GPDs. 

2.2 Probing transversely localized protons 

To discuss the spatial structure of a proton in the transverse plane we now introduce 
proton states which are transversely localized. To avoid infinities in intermediate 
steps we use wave packets 



3 $(p)b,A), (9) 

which we choose to have definite plus-momentum, $(p) = 5{p^ — p^) $_l(p). We 
find it convenient to use a Gaussian shape for $x(p), which will allow us to perform 
integrations explicitly. Abbreviating 

G(x,a^)=exp[-^] (10) 
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we thus define states centered around bo with an accuracy of order cr, 



rf2b-i^G(b-bo,a2)|p+,b,A), 



where 



b, A) = / ^ e-^^ " b+, p, A) = Imi \p+, h, A). (12) 

is completely localized in the transverse plane. Of course the proton is an extended 
object, and we will see in section IHl that more precisely it is the "transverse center 
of momentum" which is localized for \p~^, b, A) and smeared out by an amount a for 
Ip"*", b, A)o-- The normalization of these states is 

.(p'+,b',A'|p+,b,A). = p+5{p'+-p+)—^G{W-h,2a')5x'x (13) 

1 

An 

Notice that the relativistically invariant integration element in Q does not mix 
p and p^. This is not the case if the three-momentum components are taken as 
independent variables, since in (Pp dp^ / {2p^) with p'^ = (p^ + (p^)^ + m^)^/^ the 
dependence on p and p^ only decouples in limiting cases. Forming wave packets with 
definite p"*" instead of definite p^ thus keeps our expressions simple. Notice however 
that for fixed p"*" the momentum component p^ = (p"*" — p~)/\/2 becomes negative 
for very large |p| since p~ = (p^ + m?)/{2p^). If we want to have the interpretation 
of a proton moving with a well-defined longitudinal momentum, we must restrict the 
relevant range of integration in (jlll) to |p| <^p^ and thus can only have a transverse 
localization a S> l/p'*', as observed in j4j. We should thus choose a frame with large 
p'^ . This is not a restriction for our purposes, since it is just in such a frame where 
we have the natural interpretation of the proton as a bunch of partons, probed in a 
process where the GPD is measured. We also remark that if p"*" ^ m and p^ ^ |p|, 
the light-cone helicity states we are using here coincide with usual helicity states up 
to corrections of order m/p"*", see jn|. 

We now want to take the matrix elements of the operator 



0,,(z) = / ^ e"^^^- g(0, -\z-, z) 7+g(0, iz", z), (14) 



dz~ 
An 

which defines the GPDs H and between transversely localized proton states. The 
discussion for H and E is completely analogous. Using ([T]) and (fTT|) we get 



«'+ h' 



1 



p 



/tiw «(p-;^) «(p' «(^-«-^'- (15) 
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where for brevity we have omitted proton hehcity labels. Note that at this point we 
need the definition of 7i in a frame with arbitrary p and p'. Because Ti depends on 
these vectors only through D we can carry out one of the integrations after a change 
of variables, and find 



1 



1 



2\2 



e-*°-^ exp 



IGvr 



where b and 5h are given by 



167rV2 1 + ^2 



e(i-e 



D ■ 5h 



167r3 



G(5b,i(i+a^ 
2(1+ a ^ 



H(x,e,D) 
(16) 



(5b, 



b' 



6h. 



(17) 



Up to smearing by Sh, constrained by a Gaussian to be of order a, the transverse 
positions of the the initial and final state proton are not independent, and the crucial 
finding is that for nonzero ^ they are not the same. Technically this is due to the 
fact that for given nonzero ^ the matrix element Ti is not a function of p' — p but of 
D, which we saw to be a direct consequence of Lorentz invariance. More physically, 
the transverse location of the proton changes when a finite longitudinal momentum 
is transferred, which will become evident in section IHl Rewriting (fTT)|) with 5b = we 
obtain the main result of this section: 



(D b 



(27r)2 

= 647rV2 



i+e 



I _ ^2)5/2 ^ 
(1-^2)5/2 . 



p 



'+ 



p 



'+ 



1-e 



p 



1 + e 
eb 



where in the second step we have used translation invariance in the transverse plane, 
and where 

2(1+ a 1 



^ (i-e)2D2- ^^^^ 

Coming back to the freedom of choice mentioned at the end of section ITTl we remark 
that one could write p6|) to (fTHj) in terms of a different impact parameter variable, 
related to b by a .^-dependent scale factor. This would of course not change the actual 
values of the transverse distances described in these equations. 
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Figure 1: Representation of a GPD in impact parameter space. Plus-momentum 
fractions refer to the average proton momentum ^{p + p') and are indicated above or 
below lines. The region ^ < x < 1 is shown in (a), and the region \x\ < ^ in (b). 

2.3 Discussion 

We now have an interpretation of the Fourier transformed GPDs introduced in (jH)) by 
taking the limit a — of (fTHj) . Of course its right-hand side does not vanish in this 
limit, since the normalization ()13|) of our proton states includes a factor l/cr^. As data 
on GPDs will in practice only be available in a finite range of t, one might consider to 
actually perform the Fourier transform as given in p8|) . with of order of the largest 
measured value |t|max of \t\. The Gaussian damping factor would limit the effects of 
extrapolating the integrand to unmeasured values of or cutting off the integral. 
Relations Q and (fT^ give the accuracy a to which we can localize information in 
impact parameter space as a ~ (|t|max — l^lmin)"^''^ with |t|min = 4^^m^/(l — ^^). 

Figure H shows the physical picture encoded in (fTH|) . GPDs in impact parameter 
space probe partons at transverse position b, with the initial and final state proton 
localized around but shifted relative to each other by an amount of order ^b. At 
the same time, the longitudinal momenta of the protons and hadrons are specified, 
in the same way as in the t-dependent GPDs. In the DGLAP regions x G 1] and 
X G [—1,—^], the impact parameter gives the location where a quark or antiquark 
is pulled out of and put back into the proton. In the ERBL region x G [— ^,^] the 
impact parameter describes the transverse location of a quark-antiquark pair in the 
initial proton. 



7 



We remark that the shift of the transverse proton positions depends on ^ but not 
on X. The information on the transverse location of partons in the proton is therefore 
not "washed out" when GPDs are integrated over x with a weight depending on x 
and ^. As anticipated in jS] this imphes that information is already contained in 
the scattering amplitudes of hard processes, where GPDs enter through just such a 
convolution in x. Notice also that for very small ^ the difference between the proton 
positions becomes negligible compared to the impact parameter b of the struck quark. 
In this sense the situation is simpler for the impact parameter than for the longitudinal 
momentum fractions x + ^ and x — ^. Even at small their difference cannot be 
neglected when x is of order ^, which it typically is in the convolution with the hard 
scattering kernel of a physical process. 



3 The overlap representation in b-space 

A physical interpretation of momentum space GPDs has been obtained in [HI by 
writing them in terms of hadronic light-cone wave functions. We will now derive the 
analog of this in the mixed representation of definite plus-momentum and impact 
parameter. 



3.1 Fock space and wave functions 

A wave function representation of GPDs is naturally derived in the framework of 
light-cone quantization and in the gauge A'^ = 0. We will only recall essentials here; 
more detail can for instance be found in |Sj or in [Hj. In light-cone quantization the 
dynamically independent components of a quark field are given by the projection 
= ^■y~'y~^q- At a given light-cone time, which we take as = 0, this can be 
expanded in terms of creation and annihilation operators for quarks and antiquarks 
of definite momentum, 



+ d\k+, k, /i) v+{k+, fi) e'^^'~-'^-^ ] , (20) 

where /z specifies the parton light-cone helicity. Operators for quarks and antiquarks 
with definite impact parameter are given by 



r d k 

6(A;+,b,/i) = j Y^6(A;+,k,/i) e 



d{k+,h,^^) = /i!^d(fc+,k,/i)e^i'-^. (21) 
The good components of the quark field can then be rewritten as 
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[22] 



Notice that for this to work it is essential that the good spinor components m+(A;^, /i) = 
|7~7+'u(/i;, /i) and /i) = |7"7+f (A;, /i) are independent of k. In other words, for 

the good field components one can specify the light-cone helicity of a parton without 
specifying its transverse momentum. In full analogy to (|?Hl one can define anni- 
hilation operators a(A;+,b, /x) for gluons with definite impact parameter from their 
momentum space counterparts a(fc"'", k, /x). 

The Fock state expansion of a hadron state in momentum space reads 

-1 „ J™ N -1 AT AT 

ip". P. A) = E 7^ / n ^ <i - E n 'f^. (p _ ^ k, 

Nq Nq+Nq N 

X ip^p{xi,ki- Xip) Y[b^{xip'^,ki) Yl d^{xip^,ki) Y[ kj) | 0) 

i=l i=Nq + l i=Nq+Nq + l 

(23) 

using the conventions of [H]. specifies the number of partons in a Fock state, and f3 
collectively labels its parton composition and the discrete quantum numbers (flavor, 
color, helicity) for each parton. /^r/j is a normalization constant providing a factor 
n! for each subset of n partons with identical discrete quantum numbers; for ease of 
writing we have omitted labels for these quantum numbers in the creation operators. 
Notice that the wave functions depend on the momentum variables of each parton 
only via the light-cone momentum fraction Xj = kf /p^ and its transverse momentum 
kj — Xjp relative to the one of the hadron.^ 

Using the inverse transforms of (j21|) it is easy to obtain from (j23p the Fock state 
decomposition for a transversely localized proton state (jl2j) in terms of transversely 
localized partons: 

-i N J N N N 

b. A) = E ^ / n ^ ^1 - E ^0 (4^)"""' n ^'b. (b - ;^ xM. 



i=l V -^i- i=l i=l i=l 

Nq Nq+Nq N 



X ^^^{xuhi-h)Wh\xip'^,hi) \{ d^Xip"" M) H a) {xip"" ,hi) \Q) . 

i=l i=Nq+l i=N„+Nn+l 



(24) 



^Thc momentum kj — XiP is obtained from k,; by a transverse boost Q to the frame where the 
parent hadron has no transverse momentum. 
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The wave functions for definite transverse momentum or impact parameter are related 
by 



„ N N 

i,hi-h) = / [(fk]N exp[i ^ kj ■ tp^^ixi, k^) with b = ^ s^b 

i=l i=l 



i=l 

N 



ipNfsi^i, kj - Xip) = / [(fh]N exp [-i ^ kj ■ bj ip^/sixi, hi) with p = ^ kj , 

i=l 1=1 



(25) 



and normahzed to 



pA 



[dx]N [rf^kAr] ?/'^^(a;i,ki) 



[dx]N [d'^h 



N 



(26) 



where is the probabihty to find the corresponding Fock state in the proton, so 
that in total bPn0 = 1- Here we have used the shorthand notation 



[d'hU 



1 



TV 



TV 



3W- 



i=l 



i=l 



[dx] 



N 



(IGtt 

TV TV 

1=1 i=l 
TV N 

Y[dxi 6(l -^Xi 

i=l i=l 



(27) 



for the iV-parton integration elements. 

Let us briefly discuss the Fock state decomposition in impact parameter form ()24|1 . 
The parton i is transversely localized at bj, but due to translation invariance the wave 
function ip^/s depends only on its position relative to the center b of the proton. The 
constraint 

TV 

h = Y.Xihi (28) 

i=l 



identifies b as the "center of plus-momentum" PT| or "transverse center of momen- 
tum" J2] of the partons in each separate Fock state {N(3). This constraint is due to 
the invariance of the light-cone formulation under transverse boosts. Technically, it 
arises in the derivation of ()24|) from (^Hj) because the wave functions for proton states 
with different transverse momenta p have the same functional form ip^p{xi, kj — x^p). 
Notice the analogy of this situation with nonrelativistic mechanics |1H IT^. The trans- 
verse boosts (0) have the same form as a Galilean transformation in two dimensions, 
with corresponding to the mass m of a particle , and v to the velocity charac- 
terizing the transformation. The conserved quantity corresponding to the transverse 
center of momentum ()28p is the center of mass ^niiYi/ J^f^i of A^-body system. 
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3.2 The wave function representation of GPDs 



After the preparations of the preceding subsection it is easy to write the matrix 



elements Hy\ and Jyx in terms of wave functions ip^g- This can be done either by 



transforming the known representation of Tiyx in terms of V'at/?) o^' by repeating the 
derivation given in [Hj directly in the impact parameter formulation. The result in 
the DGLAP region x G 1] is 



/ 1-^ r 

Iyx{x,^,h) = EVl-e' E [dxUld'hU (29) 

N,I3 j=q 



X 6{x - X,) 5(2) (b - b,) rNp{^T\ b, - br ) i^'^Npi^. b, - bi,'^ 

with wave function arguments 



1+r ' 1-e 

in _ + ^ out _ - ^ (on\ 



and transverse locations of the proton states 



,in _ ^ 1 1 out ^ 



K = YT^b,, br = -Y3^b,- (31) 

The label j denotes the struck parton and is summed over all quarks with appropriate 
flavor in a given Fock state, and the labels {N, (3) are summed over the corresponding 
Fock states. The representation in the DGLAP region x G [—1,—^] is obtained 
from by reversing the overall sign, by changing 5{x — xj) into 6{x + xj), and by 
summing j over antiquarks. In the ERBL region x G [— ^,^] we have 

2-7V ; N 



XA'A(x,e,b)= E Vi-^ VI + ^ T.^^ (32) 

X y dxj JJ^ dxi b{\ — E ^i) 

X {Airf-' J d%, n d'h, 5(2) (eb, + E ^^b,) 

X 5{x-x,) 5(2)(b-b,) ^;v-i,,'(^r,b.-br)v^^+i,^(xr,b,-b^'^) 

with b[f and bg"* as in §^ and 

in _ ^ + in _ ^ ~ /oo\ 
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The partons j,j' are the ones emitted from the initial proton, and in (jH^ one has 
to sum over all quarks j and antiquarks j' with opposite helicities, opposite color, 
and appropriate flavor in the initial state proton, over all Fock states (A^ + 1,(3) 
containing such a qq pair, and over all Fock states (A^ — 1,(3') of the flnal state proton 
with matching quantum numbers for the spectator partons i ^ The statistical 
factors rij {riji) give the number of (anti)quarks in the Fock state (iV+ 1, (3) that have 
the same discrete quantum numbers as the (anti) quark pulled out of the target. The 
impact parameter representations for quark helicity dependent and for gluon GPDs 
are analogous to and (jH^ : the relevant modifications can easily be deduced from 
the momentum space expressions in jS]. 

If one defines the impact parameter GPDs with an additional Gaussian weight 
(7(0, D^) as in (fTHj) then one simply has to replace 

in (pnjl and (jHS))- This shows again that with data for up to order one can 
obtain information in impact parameter space to an accuracy of order 1/|Dct|. We 
remark that the precise implementation of this is different in our result ()18|). where 
the transverse location of each proton is smeared out independently, and in (jH^ . 
where the smearing is over the location of the struck parton or parton pair, with a 
smearing of the proton positions induced by (j31|) . 

The physical interpretation of the overlap formulae (j29|) and (|32p is the same as 
the one we have obtained in section 1231 and illustrated in Fig. ^ We now see why for 
^ 7^ the transverse location of the proton is different before and after the scattering: 
according to our discussion in subsection 13.11 its transverse center of momentum is 
shifted because the proton is subject to a finite transfer of plus-momentum. We thus 
find that GPDs at nonzero ^ correlate hadronic wave functions with both different 
plus-momentum fractions and different transverse positions of the partons. Notice 
however that the difference in transverse positions is a global shift in each wave 
function; the relative transverse distances between the partons in a hadron are the 
same before and after the scattering. 

3.3 Positivity constraints 

Positivity constraints for GPDs in the impact parameter representation have recently 
been considered by Pobylitsa in a very general framework JH]-^ We shall not elaborate 
on this subject in detail, but only give a set of inequalities that readily follow from 
the overlap representation (j29j) for in the DGLAP region. This representation 
has the structure of a scalar product (0^^ | 0|) in the Hilbert space of wave functions, 

■^Note that Pobylitsa defines impact parameter GPDs as Fourier transforms with respect to the 
vector (1 — 5^)D, cf. our remark at the end of section ITTl 
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with 



l-N 

l-N 



A' 

Nf3 



X 



out 



(35) 



For ^ = this imphes that X++(a;, 0, b) must be real and positive if a; > and neg- 
ative if a; < 0, as observed in ^3]. Furthermore we have the Schwartz inequahty 
1(0^^ I (f>\)\^ < (0-^ I 0-^) I which after suitable changes of integration vari- 
ables in gives 



X 



:,0 



1-r 1-e 



0, 



1+r i+e 



(36) 



for ^ < X < 1 and any combination A', A of proton helicities. Further relations are 
obtained by replacing 1 with the combinations X±X for definite parton helicity. More 
detailed inequalities involving the various proton and parton helicity combinations can 
be obtained using the matrix structure in the helicities, as shown in [T3] . 



4 Special cases 

The wave function representation makes explicit in which way GPDs probe the lon- 
gitudinal and transverse structure of a hadron in a correlated way. In particular we 
have seen that for finite ^ they are correlations of wave functions where the partons 
differ both in longitudinal momentum and in impact parameter by an amount con- 
trolled by ^. In this section we briefly discuss the type of information one obtains in 
some well-known special cases. 

4.1 The limit ^ = 

In the case ^ = (where D coincides with the transverse momentum transfer p' — p) 
the wave function representation becomes particularly simple, with 

JyA(x, 0,h) = Y.T. f [dx]N [d'hU (37) 

N,l3j=q 

X 6{x-Xj) 5(2)(b-bj) V^^^'(xi,bi)V^^^(xi,bi). 

for X > and an analogous relation for x < 0. We see that if in addition one takes 
the same spin state for the two proton states, the Fourier transform of H{x, 0, t) 
is expressed through squared wave functions and thus has a density interpretation. 
For the Fourier transform of H{x, 0, t) one obtains the difference of squared wave 
functions for a right-handed and a left-handed struck quark or antiquark. At ,^ = 
the distribution E only appears in the matrix element for proton helicity flip, and 
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its Fourier transform in X_+(x, 0,b) correlates wave functions with identical parton 
configurations but with opposite helicities of the parent proton. Giving it a probability 
interpretation is more involved ^1], and we will not elaborate on this issue. Notice 
finally that no information is obtained about E by setting ^ = in the appropriate 
matrix elements, where due to its definition (j2I) it is accompanied by at least one 
factor of ^. 

Inverting the Fourier transform (jH} we see that the usual quark density q{x) = 
H{x,0,0) can be written as the integral of X++(x, 0,b) over b, which removes the 
— hj) on the right-hand side of (jSZj)- This makes it explicit that ordinary 
parton densities specify the longitudinal momentum fraction of the struck parton, 
averaged over its transverse position in the target. 



4.2 Form factors 



Integrating H, E, H, E over x one obtains the elastic form factors of the quark 
vector and axial vector current. Since the form factors are independent of ^ one can 
evaluate their wave function representation in a frame with ^ = 0, which was the 
original choice of Drell and Yan JHl- Introducing Dirac and Pauli form factors for 
each separate quark flavor. 



ip'A'imrqiO) b,A)=M(p',A') 



Emi'+Fm 



2m 



1 /-o 

AtC j-oo 



we have impact parameter representations 
r dt jMh\^) Fl{t) 

J — oo ^ ' 



(39) 



j=q,q 

1 b'^-ib^ fO 



An Ibl 



dtJi(\hW-t 



2m 



Flit) 



= E E / [c^^bjAT 5^^) (b - bj) ilj]^'p{xi, hi) ^^^(xj, b^), 

Nfi j=q,q 

where aj is +1 for quarks and —1 for antiquarks, and A = +|, A' = — |. A cor- 
responding relation exists for the Fourier transform of the axial form factor. The 
pseudoscalar form factor decouples at ^ = and like E requires evaluation in a frame 
with nonzero ^. We see that the Fourier transform of the Dirac form factor has an 
immediate density interpretation, whereas the Pauli form factor F2 correlates wave 
functions for opposite proton helicities. To be precise, the Fourier transforms in (j39|) 
describe the transverse location of partons in a fast moving proton, irrespective of 
their longitudinal momenta 11]. This is the opposite of what we had for the usual 
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parton densities, which contain purely longitudinal information. The interpretation 
of an elastic form factor in a frame where the proton moves fast may seem unusual, 
but note that we are describing hadron structure in terms of quarks, antiquarks, and 
gluons here, so that such a frame is in fact very adequate. 

Taking higher moments in x gives form factors of operators containing deriva- 
tives. In this case, information on longitudinal structure is retained in the form of a 
momentum dependent weight. The simplest example is the quark part of the energy 
momentum tensor, whose (++) component can be written as [2] 



(p',A'|T;+(0)b,A) = uip',\' 



2m 

{p' — p)^{p' — p)^ 



m 



uip,X) (40) 



and is related with / dxxH and / dxxE by a sum rule. In a frame with ^ = the 
form factor decouples and one is left with A'^ and B^. Their representation in 
terms of momentum space wave functions has been discussed in jT7]. In impact pa- 
rameter space it is obtained from by replacing Fl with and F2 with B'^ on 
the left, and o"-,- with Xj on the right. In other words, the contributions from quarks 
and antiquarks j now come with the same sign and weighted by their momentum 
fractions. In particular, the Fourier transform of A"^ describes the transverse distri- 
bution of the longitudinal momentum carried by quarks and antiquarks of a given 
flavor. Integrating it over b one obtains the second moment / dxx{q + q) of quark 
and antiquark distributions. 



5 A tale of two scales 

So far we have suppressed the dependence of GPDs on the factorization scale fi, which 
in a physical process is provided by the hard scale Q. The dependence on this scale 
is given by well-known evolution equations [H 12 El UH] of the general form 

d f 

with evolution kernels K known up to two-loop accuracy Important for us is 

that these equations involve GPDs at the same t and ^, so that taking the Fourier 
transform with respect to D does not alter their structure. The evolution equations 
for GPDs in impact parameter space are hence the same as for their t-dependent 
counterparts. 

The physical meaning of n in our context is essentially as in the case of ordinary 
parton distributions. We recall that because of the short-distance singularities of 
QCD, light-cone wave functions and the underlying Fock state decomposition must 
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be renormalized |8j. For our purpose it is useful to understand the associated renor- 
malization scale /i as a cutoff on transverse momenta ^20^.^ Quarks and gluons in 
the presence of this cutoff are then "elementary" down to a transverse resolution of 
order — loosely speaking they have a transverse extension of that size. The wave 
functions in our overlap representations refer to partons at a given scale /i. An im- 
pact parameter GPD in the DGLAP region may thus be interpreted as describing the 
longitudinal momentum and transverse location of a quark with transverse size 
At larger values of /i, this quark may be seen as consisting of a quark plus a gluon, 
i.e., one will start to resolve its "substructure". In the ERBL region a GPD in impact 
parameter space describes a gg-pair, where quark and antiquark are each of size 
and at the same transverse position in the proton, to an accuracy again of order 

The role played by is to be contrasted with the transverse resolution a ~ 
(l^lmax — l^lmin)^^^^ discusscd in scctious 12.31 and IT^ which refers to the transverse 
location of a parton within its parent hadron. In other words, the momentum trans- 
fer t is related to where a parton is found in the proton, whereas [j? determines what 
is meant by "a parton". An ordinary quark distribution at very large /z^, for in- 
stance, contains information on the longitudinal momenta of quarks seen with very 
fine transverse resolution, but no information at all on where quarks are located in 
the transverse plane. To use an analog from optics consider a cell under a microscope. 
Then corresponds to the optical resolution and limits which details of the cell 
one can see. The analog of (|t|max — |^|min)~^'^^ specifies how precisely one controls the 
position of the cell under the objective, in order to determine where the magnified 
detail is located within the cell. 

At this point we can make a comment on the representation ()39p of the elastic form 
factors i^i(t) and F^ii) in terms of wave functions for quarks, antiquarks and gluons. 
At small values of — t, say below 1 GeV^, one may wonder how such a description 
is possible, given that the form factors are measured in elastic scattering processes 
whose momentum transfer t is insufficient to resolve any partonic structure at all. 
The solution of this apparent paradox is that these form factors are independent of 
a renormalization scale //^ because the quark vector current (l{z)'y^q{z) is conserved. 
They are hence the same when evaluated at /i^ ~ —t or at /z^ of several GeV^. Phys- 
ically speaking, the transverse distribution of charge in the proton is independent of 
how much substructure is resolved, so that one may represent this charge distribution 
as due to quarks and antiquarks, even though one has not resolved them explicitly. 
In this sense, Fi{t) and F2{t) at small t do contain information on partons, but this 
information is not specific to these degrees of freedom. 

Notice that the situation is different for form factors of other operators, which 
are given by higher moments of GPDs in x. The quark energy-momentum tensor 
for instance does depend on the renormalization scale, and the information of its 

''Of course such a cutoff regularization — which among other things breaks Lorentz invariance — 
has its limitations, but it does make the physics transparent. There are other possible cutoff schemes, 
involving for instance light-cone energy fc~ instead of transverse momentum 
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form factors is specific to the value of /x. As n increases, the average momentum 
carried by quarks at a given impact parameter will become smaller, since one resolves 
processes where the quarks lose momentum by radiating gluons. In the case of energy- 
momentum one still has a conserved current when summing over all quark flavors 
and gluons, but this does not hold for other local operators connected to moments of 
GPDs. 

6 Unintegrated parton distributions 

Apart from GPDs, there is another class of nonperturbative functions that carry 
information not only on longitudinal but also on transverse hadron structure. These 
are fcr-dependent or unintegrated parton distributions. Let us see how the information 
they contain looks like when expressed in terms of transverse position, and contrast it 
with the picture we have obtained for GPDs. For simplicity we will restrict ourselves 
to forward distributions and set p = p'. It is then sufficient to consider a frame with 
p = 0. Following the naming scheme of |23 ^6 define 

k) = I e.xp+.--.k.z ^1 -(0^ ^+^(0^ 1^-^ (42) 

which is related to the usual quark density by / (Pk k) = q{x). We suppress again 
the dependence on the factorization scale /x, whose discussion proceeds in analogy with 
the preceding section. A word of warning is in order concerning the absence of the 
usual Wilson line between the operators q and q. These are now taken at different 
transverse positions ±|z, so that a Wilson line appears even in the gauge A'^ = 0, 
except for particular choices of the integration path. As recently discussed in [22] there 
is subtle physics encoded in the choice of path and the Wilson line. This is beyond 
the purpose of our present investigation, which is a basic understanding of the spacial 
information contained in various types of parton distributions. The same holds for 
recent work where it was argued that even for the usual parton distributions, with 
z = 0, physical effects of the Wilson line are not correctly reproduced in light-cone 
gauge pi. 

In transverse momentum space, the distribution /i(a;, k) can be interpreted as a 
probability density. This is readily seen from its wave function representation, which 
involves kj)p, with the wave function arguments of the struck parton fixed to 

be Xj = X and = k.^ In terms of impact parameter wave functions we have 

/ d^k e^i'-^ k) = ^ ^ / [dxU [d'^hU (43) 

X S{X- X,) rN(3{x., br - br ) ^U^^, K - K) 

^In a frame with p 7^ one has instead kj = k — xp, in accordance with the invariance of kj 
under transverse boosts ©. 
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Figure 2: Impact parameter representation of an unintegrated parton distribution, 
z is the Fourier conjugate variable to the transverse momentum k of the struck parton. 



for X > 0, with wave function arguments 

= br = b, fbr z ^ J, 

= b, + ^, b°- = b, - ^, (44) 
and proton states centered at 

K = br = -^l- (45) 

The corresponding physical picture is represented in Fig. |2l As is already evident 
from the operator defining k) in (jl^ . we see that the struck quark is at a 
relative transverse distance z in the initial and the final state. Unintegrated parton 
distributions thus describe the correlation in transverse position of a single parton. 
Notice that, in contrast to the situation for GPDs, the struck quark now has a different 
transverse location relative to the spectator partons in the initial and the final state 
wave functions, in addition to the overall shift of the proton center of momentum 
dictated again by Lorentz invariance. The dependence of the Fourier transform (j43|) 
on z thus describes how much the transverse location of a single parton can vary in 
the proton wave function when all other partons are kept fixed. Notice finally that 
since /i(x,k) is measured at t = 0, the overall position of the struck quark with 
respect to the center of the proton is integrated over in ()4Hj) . 

The quantities that contain this information in addition are of course /ct unin- 
tegrated GPDs. They have been used in an estimate of power corrections for hard 
exclusive processes in ^| and recently been discussed in the context of small-x physics 
in [23] . Their impact parameter representation combines the characteristic features of 
fl29|) and (jH^ with those of ()43|). In the DGLAP region, both the average transverse 
position b and the relative offset z of the struck parton in the initial and final state 
are now specified, whereas the ERBL region describes a gg-pair in the initial state 
proton, with an average transverse position b and a relative separation z between 
quark and antiquark. 
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7 Conclusions 



A natural setting to represent the structure of a hadron in terms of quarks and gluons 
is a frame where the hadron moves fast. Its direction of motion and the two perpen- 
dicular ones then acquire different physical roles, and a comprehensive description of 
the hadron requires information both on the longitudinal and the transverse degrees 
of freedom. Generalized parton distributions are among the quantities that allow one 
to address such questions as: "How broad is the spatial transverse distribution of fast 
quarks compared with slow ones, and how does it compare with the spatial transverse 
distribution of gluons?" 

Whereas the description of scattering processes is in momentum space, useful phys- 
ical insight can be obtained in the mixed representation of longitudinal momentum 
and transverse position, or impact parameter. This representation is natural in var- 
ious contexts, for instance in the description of high-energy scattering processes, see 
e.g. |2B], or the resummation of Sudakov logarithms in hard reactions Burkardt 
has pointed out that for zero skewness ^, GPDs in impact parameter space have the 
simple representation of a joint density in the longitudinal momentum of a parton and 
its transverse distribution in the proton, at least for helicity conserving transitions 
|3]. We find that the situation for nonzero ^ is similar up to an important difference: 
as the proton loses longitudinal momentum its transverse position is shifted by an 
amount proportional to ^. This is because as a consequence of Lorentz invariance 
the "transverse position" of the proton is the vector sum of the transverse positions 
of its partons weighted by their longitudinal momentum fractions. A more detailed 
picture is obtained by expressing GPDs in terms of impact parameter wave functions, 
related to the ones in momentum space by a Fourier transform. In addition to their 
longitudinal momentum fractions, the transverse distance of partons from the proton 
center differs in the initial and final state, but their relative distance to each other in a 
hadron stays the same. This is in contrast to the case of fcr-dependent (but forward) 
parton distributions, which in impact parameter space are correlation functions for 
the transverse distance of a single parton with respect to all other partons in the wave 
function. 

All together, GPDs at nonzero ^ describe quantum mechanical correlations within 
a hadron rather than probability densities, both in the transverse momentum and the 
impact parameter representations. Despite their complexity in detail, the underlying 
physical picture is simple, as shown in Fig. ^ As to achievable spatial resolution, 
GPDs measured in a process with hard scale provides information on partons as 
seen with resolution Q~^, whereas the invariant momentum transfer \t\ to the proton 
must be known from its minimum value |t|min up to |t|max in order to locate partons 
in the transverse plane to an accuracy of order (|t|max — |^|mm)~^^^- 
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Erratum 



Generalized parton distributions in impact parameter space 

Eur. Phys. J. C25, 223 (2002), |hep^h/0205208 
M. Diehl 



The kinematical factor in the positivity bound (36) is incorrect. The bound cor- 
rectly reads 

(1 - er\iMx,^M f < x++(^,o, ^) x++(^,o, ^). 

Our corrected result agrees with inequality (25) in |13j, taking into account the dif- 
ferent normalization conventions here and there. 
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